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Abstract 

The thermodynamic Bethe ansatz (TBA) and the excited state 
TBA equations for an integrable spin chain related to the Lie superal- 
gebra osp(l\2) are proposed by the quantum transfer matrix (QTM) 
method. We introduce the fusion hierarchy of the QTM and derive the 
functional relations among them (T-system) and their certain combi- 
nations (K-system). Their analytical property leads to the non- linear 
integral equations which describe the free energy and the correlation 
length at any finite temperatures. With regard to the free energy, 
they coincide with the TBA equation based on the string hypothesis. 
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1 Introduction 

Solvable lattice models related to the Lie superalgebras |jj have attracted 
a great deal of attention [|, |3], [|, ||, H S S 3- For example, the super- 
symmetric t — J model in strongly correlated electron system has received 
much attentions in relation with the high T c superconductivity. These mod- 
els have both fermionic and bosonic degree of freedom, and are given as 
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solutions of the graded Yang-Baxter equations 0. To solve such mod- 
els, the Bethe ansatz is widely used in many literatures (see, for example, 
0, 0, |TT], 0, |13], |IJ, |IJ, |I§ 0, |IJ, [19], [2(J and references therein). However, 
many of them deal only with models related to simple representations like 
fundamental ones; there was few systematic study by the Bethe ansatz on 
more complicated models such as fusion models ||21|| . 



In view of such situations, we have recently carried out [^, £3], |24|, |25|, £| 
systematically an analytic Bethe ansatz [^7], [ft], [31], [ftj, |33| related 

to the Lie superalgebras sl(r + l\s + 1), B(r\s), C(s) and D(r\s). Namely, 
we have proposed a class of dressed vacuum forms (DVFs) labeled by Young 
(super) diagrams and a set of fusion relations (T-system) among them. 

Besides the eigenvalue formulae of the transfer matrices, the thermody- 
namics have also been discussed by several people. Particularly, the thermo- 
dynamic Bethe ansatz (TBA) []33| equation was proposed for the supersym- 
metric t — J model fllQl , which is related to s/(l|2), and the supersymmetric 
extended Hubbard model Jill, which is related to s/(2|2). Moreover, there 
is a paper [35] on the excited state TBA equation for these two important 
models of ID highly correlated electron systems from the point of view of the 
quantum transfer matrix (QTM) method. In addition, the TBA equation for 
sl(r\s) model was presented J36J in relation with the continuum limit of the 
integrable super spin chains. 

However, the thermodynamics of the quantum spin model related to the 
orthosymplectic Lie superalgebra osp(r\2s) is not so understood as sl(r\s) 
case. In particular, as far as we know, there have been no literatures on the 
TBA equation even for the simplest orthosymplectic osp(l\2) integrable spin 
chain 0, f|, |^, [7|, |8|, [|, [I7|| . This is regrettable because this model may 
be related to interesting topics such as N = 1 superconformal-symmetry in 
field theory, and the loop model, which will describe statistical properties of 
polymers in condensed matter physics In view of these situation, we 

have recently proposed the TBA equation for the osp(l\2) integrable spin 



chain |38] by using the string hypothesis |39], [40], f4~T| 



Though we expect that the resultant TBA equation describes the free 
energy correctly, there exists the deviation from the string hypothesis 



42, 



43 



In addition, it is difficult to evaluate other physical quantities such 
as correlation length. 

As an alternative method which overcomes such difficulties, the QTM 
method has been proposed ||6], [48], [4^, [50], [5T| . Now we shall briefly 



sketch the QTM method. Utilizing the general equivalence theorem [|45|, [46 



one transforms the ID quantum system into the 2D classical counterpart 
and defines the QTM on such a fictitious system of size N (referred to as the 
Trotter number, which should be taken iV — > oo). Since the QTM has a finite 
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gap, the original problem for the calculation of the partition function reduces 
to finding the single largest eigenvalue of the QTM. To evaluate it actually, 
we utilize the underlying integrable structure, which admits introduction of 



the "commuting" QTM with a complex parameter v [p2| , [53], [54], |5^, [56], [57| . 
Furthermore, we introduce some auxiliary functions including the QTM itself, 
which satisfy functional relations. 

We select these auxiliary functions such that they are Analytic, iVonZero 
and have Constant asymptotics in appropriate strips on the complex f-plane 
(we call this property ANZC). Thus we can transform the functional relations 
into the non-linear integral equations (NLIE) which describe the free energy. 
In these NLIE, we can take the Trotter limit N — > oo. 

Adopting a subset of fusion hierarchy as auxiliary functions, we find that 
these NLIE are equivalent to the TBA equation based on the string hy- 



pothesis [^, |2|, 0, In general, a set of fusion hierarchy satisfies 

the functional relations called T-system, and this T-system is transformed 
to the V-system |. By selecting them so that they have ANZC property 
in appropriate strips, one derives the NLIE which will be identical to the 
TBA equation. Furthermore, considering the sub-leading eigenvalues of the 
QTM, we can derive systematically the "excited state" TBA equations which 
provide the correlation length at any finite temperatures. 

The purpose of this paper is to apply our recent results PBf to the osp(l\2) 



integrable spin chain, and to construct the TBA equation and it's excited 
state version from the point of view of the above-mentioned QTM method. 
We have also confirmed the fact that our TBA equation coincides with the 
one |38| from the string hypothesis. We believe that this paper yields a 
basis of future studies of the thermodynamics by the QTM method for more 
general models such as the osp(r\2s) model. 

The layout of this paper is as follows. In section 2 we formulate the 
osp(l\2) integrable spin chain at finite temperatures in terms of the com- 
muting QTM. In section 3, to evaluate its eigenvalue, we utilize the fusion 
hierarchy of the QTM. It is given as a set of the dressed vacuum forms (DVFs) 
which is a summation over tableaux labeled by a Young (super) diagram with 
one column. These DVFs and their certain combinations (V-functions) sat- 
isfy the T-system and the F-system, respectively. The formulation of the 
DVFs and the functional relations are essentially independent of the vacuum 
part. Thus we can utilize the results in Ref. p6[ only by replacing the vac- 



uum parts of the DVFs defined in Ref. [Ell with those of the QTM. Based 



on the analytical property of the DVFs and the V-functions, in section 4 we 



a For simple Lie algebras, general relations between T and F-system are given in 
and || (see also @). 
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derive the non-linear integral equations for the free energy, and in section 
5 for the correlation length. Section 6 is devoted to the summary and the 
discussion. 



2 Quantum transfer matrix 

In this section, we introduce the quantum spin chain related to the Lie su- 
peralgebra osp(l\2) and formulate the commuting QTM for this model. 

The classical counterpart of this model is given by the ^-matrix (see for 
example, Ref. ||) 

R(v) = I + vP 9 —^-E, (2.1) 

which is a rational solution of the graded Yang-Baxter equation [Q] associated 
with the three dimensional representation of the Lie superalgebra osp(l\2) 
(to be precise, "super Yangian" Y(osp(l\2))), whose basis is Z 2 -graded and 
labeled by the parameters p(l) = p(T) = 1, p(0) = 0. Here P 9 denotes 
the graded permutation operator which has 9x9 matrix elements (P 9 ) a f = 
(—iy( a M b )fi a d fi b c anc i e j s a g x 9 matrix whose matrix elements are = 
Ci.abi.Oi~ 1 ) cd where a is a 3 x 3 matrix 

a n «i an 

a = | ct i «oo "oT | = | 1 | , (2.2) 
"Ti "To "TT 

and {a,b,c,d} G {1,0,1} with the total order 1 -< -< 1. The row-to-row 
transfer matrix T(v) is defined by 

T(v) = Ti a [R aL (v)---R al (v)}, (2.3) 

where L is the number of lattice sites and R a j (v) denotes 

R(v) = PR(v), (2.4) 

which acts non-trivially on the auxiliary space a and the j-th site of the 
quantum space. Note that P°£ = 5 a ^b, c is the (non-graded) permutation 
operator. The Hamiltonian of the corresponding quantum system is given 
by taking the logarithmic derivative of above transfer matrix (|2.3|) at v — 0, 





f o 







H 





i 


i 




V-i 








7 = 1 N ' 



(2.5) 
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where we assume a periodic boundary condition. Here J is a real coupling 
constant which determines the phase of this model; the ferromagnetic and 
antiferromagnetic regimes correspond J > and J < 0, respectively (see for 
example, |17|, |38|). 

To consider the finite temperature property of the model (|2.5| ), we shall 
introduce another transfer matrix T(v) constructed by the i?-matrix R(v), 
which is defined by 90° rotation of R(v), i.e., Rjkiv) = tk Rkj[v) (£& means 
the transposition for the i?-matrix in the fc-th space): 

f(v) = Ti a [R aL (v)---R al (v)}. (2.6) 



We can see that the logarithmic derivative of above transfer matrix ( |2.6| ) 
also represents the Hamiltonian fl2.5|) . Thus the expansion of the transfer 
matrices ( |2.$ ) and ( |2.6| ) are expressed as 



T(v) 



T(0)|l + ^ + O( W 2 )| 



T(v) = T(0) jl + yt; + eV)j. (2.7) 

Combining above relations and using the fact T(0)T(0) = 1 (note that T(0) 
and T(0) denote the right and left sift operators, respectively), we have 

T(v)f(v) = l + ^v + 0(v 2 ). (2.8) 



Consequently, the partition function Z of the model (|2.5|) can be expressed 



as 



Z = Tre"^ = lim Tr (t(u n )T(u n )) N/2 , u N = (2.9) 

where N is an even number: it is called the Trotter number, which represents 
the number of fictitious sites in the Trotter direction; f3 = l/(fcsT) (T is the 
temperature). Thus the free energy per site / is given by 

1 ~ ~ 

/ = - lim lim y^lnTr (t(u n )T(u n )) 2 . (2.10) 

However, in the antiferromagnetic case J < 0, the eigenvalues of T{un)T{un) 
will be infinitely degenerate in the limit N — ► oo. Therefore taking the trace 
in above expression (|2.10|) is a serious problem. To avoid this difficulty, we 
transform the term T{un)T{un) as follows: 

N N/2 

Tr (t{u n )T{u n )Y = TrJ] 

a 2k, a 2k-l 

k=l 
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L N/2 

= TrYlTi^YlRa^jiu^Ra^jiuN). (2.11) 

j=l k=l 

Now we introduce the QTM, which plays the fundamental role to describe 
the thermal quantities 



TP(u,v) = TTj Y[ Ra 2m j{u - iv)R a2m _ ld (u + 



N/2 



IV 



(2.12) 



m=l 



Due to the Yang-Baxter equation, we see that the QTM is commutative: 

{Ti l \u,v),Ti l \u,v')}=0. (2.13) 

Hereafter we write the fc-th largest eigenvalue of the QTM T^(u,0) as 
Tj"2(w, 0). Since the two limits in ( |2.10| ) will be exchangeable as proved in 
|f46| , |47|| , we shall take the limit L — > oo first. Since there is a finite gap 
between t[i(un, 0) and (un, 0) even in the Trotter limit N — > oo, we 
have 

(2.14) 



/ = -- lim lnT^^O). 

jj N-*oo 



The thermodynamical completeness — lim^o Pf — hi 3 follows from (0, 0) 
3, which is clear from #12(0) = P12. 

Taking the ratio of the largest eigenvalue and the second largest one, we 
can systematically calculate the correlation length £: 



- = — lim In 



Tfi,W,0) 



7uW,0) 



(2.15) 



3 Fusion hierarchy and functional relations 

In this section, we present the eigenvalue formulae of the fusion QTM (u, v) 
and the functional relations among them by using the results for the row- 
to-row transfer matrix |[26|| . Although the quantity we want to evaluate is 
only T±(u, 0), consideration on all {Tm\u, v)} plays an essential role in our 

b We also use the expression Tm (u, v) as the eigenvalues of the fusion QTM. 
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formulation. For a G {1, 0, 1}, we define the function [a\ i as 

Q(v-I) 



111 = Mv) 



Q{v+iv 

Q{y)Q{y + z i) 
Q{v + i)Q{v + i) 
Q(v + 2i) 
Qiy + i) ' 



(3.1) 



where Q(v) = YYj=i( v ~ v j)'i w 6 {0, 1, ... , N} is a quantum number; Vj G C. 
The vacuum parts of the functions [a\ , ( |3.1| ) are given as follows: 



^i(u) 
^o(«) 



1) 



AT— n 



V- 



(«"?) 



</>+(v)0_(», 

, N - n <t>-( v )<t>+(v - + £ 

1 j «M* + f) 



(3.2) 



where 



0±(f ) = (u ± mja". 



(3.3) 



The eigenvalue formula of the QTM is given as the summation over these 
boxes (O): 



(3.4) 



We have f] 

# es „ = _^([^+[T]j = o, 

under the following type of the Bethe ansatz equation (BAE): 

<t>-(vk + %)(/>+(vk -i) _ _,_ 1 , N - n Q(vk ~ %)Q(vk + i) 



(3.5) 



<P-{vk - h)<f>+(v k - 2i) 



Q(v k + ~)Q{v k - i) 
for k G {1,2, . . . ,n}. 



(3.6) 



Here Res v=a f(v) denotes the residue of a function f(v) at v = a. 
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Thus (u,v) (pi4| ) is free of polesQ under this BAE (|3.6|) . We observe, from 
numerical analysis, that the largest and the second largest eigenvalues lie in 
the sector n = N and n = N — 1, respectively. 

We shall introduce a function Tm (u,v) with a spectral parameter v E C 
and a Young (super) diagram [] (l m ), which is a candidate of the DVF for the 
fusion QTM. From now on, we often abbreviate the parameter u to simplify 



the notation. For m E Z>x, 
of the boxes in (13.11) : 



we define Tm\v) as a summation over products 



AC 



{i fc }GS(l m ) 



1 2 



-1 + m-3 



(3.7) 



-m+l 



where the spectral parameter t> is shifted as t> + . . 

'"'^ is a set of tableaux (z& e {1,0, 1}) 



from the top to the bottom. B{Y 
obeying the following rule (admissibility conditions) 



2i 



h ■< «2 d 



(3i 



The normalization function N m (v) is introduced so that the degree of T^(v) 
is 2iV with respect to v. 



AC 



TO — 1 

n 



m + l — 2j 
v H i 



m + 1 — 2j 



(3.9) 



The DVF ( 
that t2 ] (v 



is free of poles under the BAE ( |3.6|) for any m E Z>x- Note 
for m = 1 reduces to (|3.4j ) since J\fi(v) = 1. The BAE 
( |3.6| ) for the QTM has the same form as the one ^ |TB|, ^] for the row 
to row transfer matrix except the vacuum part. Thus the dress part of the 



mi.AiH.il j.~ro^(t>) ( |3.7| ) is same as that in PB| . Therefore the DVFs Tm\vj 
(m E Z>i) satisfy the following functional relations [| (osp(l|2) version of the 
fusion relations or the T-system | ), which has essentially the same form in 



d Here singularities of the vacuum parts of the DVFs, which can be removed by multi- 
plying overall vacuum functions are out of the question. 

e Note that this diagram corresponds to the Kac-Dynkin label 2m. 

f Of course, there is a large class of functional relations among DVFs labeled by any 
skew- Young diagrams. The functional relations, which we are now dealing with is nothing 
but a subset of it. However we know, from our experience, that the functional relations 
among DVFs labeled by rectangular Young diagram often have a good analyticity relevant 
for the TBA. The T-system fl3T0| ) is one of them. 

s In this paper, we use the variable m in (3.1C) slightly different from that in p6| for 
practical reason: the variable m in fl3.10| ) corresponds to 2m in |26fl . 
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]2§: 

Iff (* - ^ ) TW ( t; + I ) = eWlfiiW + TgHv)TW(v), (3.10) 
where 

(v + ^i)0 + O - - f )0+(v + f ) 



for m e Z>i. (3.1i; 



We can prove these functional relations by the Jacobi identity and a duality 
among the DVFs [p6[| . 

Using Tm (v) (P-7|), we introduce the following functions ( V-functions) : 



Y m (v):= ?i| for ( 3 - 12 ) 



Tm{v)Tm{v) 

These functions ( |3.12 ) satisfy a so-called V-system: 



lU«--)y m (t; + -) = i + { y m(v)} -i for ^eZ> l5 (3.13) 

where T_l(v) = Y (v) = 0. In the subsequent sections, these functional 
relations (|3.13|) will be transformed into the TBA equations. 



4 Thermodynamic Bethe ansatz equation 

Let {T^\(u, v)} be the DVFs {Tm(u, v)} constructed by the roots of BAE 
( |3.6| ) which provide the fe-th largest eigenvalue of the QTM T± (u, v ), and 
{Y mi k(u, v)} be the V-functions constructed by {T^\(u,v)} as in ( |3.12| ). 

In this section, we study the analyticity of the DVFs {T^\(u,v)} and 
the V-functions {Y m ^(u,v)} in the complex v -plane. Then we derive the 
non-linear integral equations (NLIE) which provide the free energy. For this 
purpose, by keeping the Trotter number N finite, we have performed numer- 
ical analyses with various values of (3 and N in investigating the location of 
zeros of {T^\(u,v)}. As seen in Fig.|], there are N/2 two-string solutions of 
the BAE (|3.6|) , which will provide the largest eigenvalue of the QTM. These 
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roots locate symmetrically with respect to the imaginary axis and the line 
%v = 3/4. 

From the definition of the DVFs ( |3.7| ), they have 2N zeros on the complex 
v-plane. For example, the location of zeros of T^\(u, v) for u = 0.05, N = 12 

is plotted in Fig.||. According to them, we observe that T^\(u, v ) have N 
zeros on the smooth curve near the line = ±(m + l)/2 and N zeros on 
the smooth curve near the line = ±(m/2 + 1). Thus we expect that the 
following conjecture is valid even in the limit iV — > oo. 

Conjecture 4.1 All the zeros of T^\(u,v) for small \u\ <C 1 are located 
outside of the strip G [—1/2, 1/2]. 



From the definition (|3.7| ), one observes the fact that all the zeros are sym- 
metric with respect to both real and imaginary axes. As seen in the figure, 
the deviation from the line is very small, which will be smaller as u — > 0. 

Once this conjecture is assumed, we can identify the strip where the Y- 
functions {Y ml (u,v)} have the property of Analytic iVonZero and Constant 
asymptotics in the limit v — > oo (we call it ANZC property). Using the 
definition of the DVFs ( |3.7| ) and the F-functions (|3.12| ), we can obtain the 
asymptotic value of the ^-functions in the limit \v\ — > oo: 

l im YmAv) = ^L±A. (4.1) 

\v\—>oo 2 

From conjecture P~T|and above asymptotics ( |4.1| ), the functions 1 + Y m> i(u,v) 
and 1 + {Y mt i(u, have ANZC property in the strip G [—£,£] (0 < 

e< 1). On the other hand, the functions Y m i(u, v) have the ANZC property 
in the strip G [—1/2,1/2] (we call this strip physical strip), except for 
Yi t x(u,v) possesses poles (resp. zeros) of order N/2 at ±(1/2 + u)i (resp. 
±(1/2 — u)i) in the physical strip for J > (resp. J < 0). Note that u = 
is a small quantity given in fl2.9|) . Using Cauchy's theorem, we can transform 
the F-system ( |3.13j ) into the NLIE in the following way. Let us consider the 



Y-system ( j3.13j ) for m > 2. First, we take the logarithmic derivative and 



perform the Fourier transformation on both side of the Eq. ( |3.13| ). Second, 



by Cauchy's theorem, the Fourier integral for the logarithmic derivative of 
Y m ,i{v) is represented by that of 1 + Y m±11 (t>) and 1 + {F ml (t>)} _1 . Fi- 
nally, performing the inverse Fourier transformation and integrating over v, 
we obtain the desired NLIE. The integral constants are determined by the 
asymptotic value in ( |4.1| ). In the case m = 1 in ( |3.13| ), we have to modify 
Yi t i(v) to transform F-system into the NLIE since Yi^(v) have the poles (or 
zeros) in the physical strip. 

Y Xjl {v) = Y hl (v){t&nh^(v + t(l/2±u))t&nh^(v -i(l/2±u))} ±N/2 , (4.2) 
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where the + and — signs in front of u and N/2 should be chosen according 
as J > and J < 0, respectively. According to the identity tanh|(w + 
i) tanh | (i> — i) = 1, li,i(i> ±i/2) in the lhs of the Eq.( ^.13[ ) can be replaced 
by Y\ t i(y) in (|4.2|) . Thus the F-system ( p. 13 ) can be transformed to the 



NLIE in the similar way as mentioned above. Consequently the resultant 
NLIE are represented as 

N ( IT TT ~) 

hi ^(v) = =Fy ln|tanh-(f + i(l/2±u))tanh-(f -i(l/2±w))| 

+K * ln(l + Y 2>1 )(v) - K * ln(l + Y^)(v), 
lnY m>1 {v) = K*hi(l + Y m+ltl ){v) + K*lji(l + Y m ^ 1 )(v) 

-K*\n{l + Y-\){v) formGZ> 2 , (4.3) 

where A * B(v) denotes the convolution 

/oo 
A(v - v')B(v')dv', (4.4) 
-00 



and 



K(v) = : . (4.5) 

v ; 2 cosh 7™ v ; 



In above NLIE ( |4.3j ), the Trotter limit N —> 00 can be calculated: 

N ( 7T 7T 

=F lim — In <^ tanh — (v + i(l/2 ± u N )) tanh — (v — z(l/2 ± u N )) 
iv— >oo 2 12 2 

TT/3J 



COsh 7TV 



(4.6) 



We thus arrive at the NLIE for Y m> i (v) which are independent of the Trotter 
number N. 

lnYii^) = *P J +K* ln(l + Y 2>1 )(v) - K * ln(l + Y^ 1 ) (v) , 

COSn 7TV 

lny m> i(r;) = K * ln(l + r m+1 ,i) (v) + K * ln(l + r m -i,i)W 

-ir*ln(l + l^i)(t;) formGZ> 2 , (4.7) 



Under the identification 1^,1 (f) = r) m (v), above equations (|4.7|) are identical 
to the TBA equation (30) in P5| , which was derived by the thermodynamic 
Bethe ansatz based on the string hypothesis. From above equations and the 
asymptotics (|4.1|), we can determine the F-functions Y m>1 (v) uniquely. 
To obtain the free energy per site, we shall modify T^(v) as 
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From the definition of the F-function (|3.12 ) and the T-system fl3.10| ) T{ 
satisfies the following functional relation. 



f$(v) 



F(»)(l + l'uW), 



where 



F(v) 



<f> + (v + i)d>.(v - |) 



ANZC property of the both sides leads to 

]nT${u N ,v) = G*]n(l + Y ltl )(v)+\n<l) + (v-i) 



(4.9) 



(4.10) 



(v + i) 



+N 

where G(v) denotes the kernel, 

G(v) 



2e 2 sinh kujy cos kv 



k(2 cosh 



£-1) 



dk. 



(4.11) 



2 sinh |7Tf 
sinh 2irv 



(4.12) 



Calculating the Trotter limit iV — > oo, we obtain the free energy per site 



47T 



\3V3 



1 ) — ksT / ln(l + Y 1A (v))dv. 



(4.13) 



This representation coincides with Eq.(33) in p8 



5 Excited state TBA equation 

In this section, we develop the analysis for the largest eigenvalue of the 
QTM to the the analysis of the second largest one. Namely we derive the 
NLIE (excited state TBA equation) which provide the correlation length at 
finite temperature. We assume that the coupling constant J is negative: 
we consider only the antiferromagnetic regime in this section. After some 
numerical analyses, we observe that n = N — 1 roots {v^} of the BAE ( |3.6|) , 
which will bring us the second largest eigenvalue of the QTM, form y — 2 
two-strings and one three-string if iV e 4Z (cf. Fig. ||); — 1 two-strings and 
one one-string if iV G 4Z + 2 (cf. Fig. |5|). As in the largest eigenvalue case, 
the distribution of these roots are symmetric with respect to the imaginary 
axis and the line $sv = 3/4. 
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From the definition of the DVF 



has 2N zeros on the 
(i) 



complex f-plain. From the numerical analyses, we observe that T^ 2 [u, v 
has N — 2 zeros on the smooth curve near the line ^sv = ±(m + l)/2, N — 2 
zeros on the smooth curve near the line = ±(m/2 + 1), two zeros on 
the imaginary axis near the points v = ±(2m + 3)i/4 and two zeros ±x m 



> 0) on the real axis. For example, the location of zeros of T 7 



(i) 



for u = 0.05 and N = 12, 14 are plotted in Figs. f||| We expect that the 
following conjecture is valid even in the limit N — > oo. 



Conjecture 5.1 (u,v) has two real zeros ±x m . 
located outside of the physical strip Qx> G [—1/2, 1/2]. 



Every other zero is 



As in the previous section, we shall identify the strip where the Y- 
functions {Y mt2 (u,v)} have the ANZC property. From the definition of 
the DVFs (| 



of Y mt2 (v) 



and the F-functions (|3.12|) , we find the asymptotic value 

[2m + 3) - 3 



\v\ 



lim Y m>2 (v) 



1 



4 



(s.i; 



From the definition of (|3.12|) and conjecture |5TT| , we find (1 + F m+ i )2 (t>))(l + 
Y m ^ 2 {v))/{l + {Y m ^{v)}- 1 ) has the ANZC property in the strip G [-e,e] 
(0 < e <C 1), while Y m ^(v) has zeros of order 1 at v — ±x m+ x, ±x m -i and 
poles of order 1 at v — ±x m if m G Z> 2 ; Y\p(v) has zeros of order N/2 
at v — ±(1/2 — u)i, zeros of order 1 at v — ±x 2 and poles of order 1 at 
v = ixx, and that every other singular point locates outside of the physical 
strip. Thus we have to modify Y m)2 (v) as follows: 



Y m o(v) 



f 7T 1 7T 1 

Y m ,2(v) < tanh -(v + i(~ - u)) tanh -(v - i(- - u)) 

x jtanh^O - x m ) tanh^(w + x m )| 

x | tanh - (v - x m+1 ) tanh - (t> + a; TO+1 ) J 



°m,l 



x <j tanh ^(v - z m _i) tanh ^(w + x m _i] 



-l + Sr, 



(5.2) 



These Y m ^ 2 (v) ( |5.2j ) have the ANZC property in the physical strip. After 
similar calculation in the previous section, we obtain the excited state TBA 
equation: 



lnYi» 



cosh 71V 



+ K * In 



1 + F 2 



2.2 



1 + ^72 
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( 7T 7T - 

— In jtanh — (t> — sci) tanh — (t> + si 

f 7T 7T 

+ In jtanh — (t> — x^) tanh — (u + X2 



\nY„ 



K * In 



v l + 5^m+l,2j (1 + Yrn-1,2 



y-1 
1 m,2 




+ 717, 



— In jtanh — (v — x m ) tanh — (v + x m ) j 
+ In jtanh -(w - x m+ i) tanh -(t> + x m+ i / 



7T 



+ In <; tanh — (u 



7T 



.i)tanh-(u + x m _i ; 
for m e Z> 2 , 



+ 



1 + 



-Til 

(5.3) 



where the integral kernel is defined in ( |4.5|) . In addition x m is determined 
by the condition Y m {x m =ti/2) = —1. Through the above excited state TBA 
equation (|5.3j), the second largest eigenvalue is given as follows: 

/ 4tt 



lim In 



-PJ 



V3\/3 



1J + 2 In (tanh ^ 



+ 



G{v) ln{(l + Yifl(v)) tanh ^(v - sx) tanh ^(v + xi)}dv.(5A) 



This second largest eigenvalue ( |5.4j ) together with the largest one ( |4.13 ) de- 
scribe the correlation length at any finite temperature (see, (|2.15|) ). These 



excited state TBA equation ( |5.3j ) and Eq. ( |5.4| ) will be difficult to obtain by 
the string hypothesis; they exemplify the efficiency of the QTM method. 



6 Summary and discussion 

We have applied the QTM method to the osp(l\2) integrable spin chain. 
Making use of the T-system and the F-system, we derived the infinite set of 
NLIE which provide the free energy and the correlation length at any finite 
temperatures. As concerns the free energy, these NLIE are identical to the 
TBA equation based on the string hypothesis. One can say that the validity 
of the osp(l\2) T-system in Ref. |26| has been confirmed from the point of 
view of the QTM method. To the author's knowledge, this paper is the first 
explicit derivation of the excited state TBA equation for the osp(l\2) spin 
chain. 

We comment on the ferromagnetic regime J > in the second largest 
eigenvalue sector, which we have not considered in section 5. In the an- 
tiferromagnetic regime J < 0, the correlation length are characterized by 
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the additional zeros of the fusion QTMs. These zeros are real and symmetric 
with respect to the imaginary axis and this symmetry will never break at any 
finite temperature. While in the ferromagnetic regime (J > 0), the "level 
crossing" will occur successively This attributes to the change of distribution 
patterns of the additional zeros as follows. First at high temperature, two 
additional zeros are on the imaginary axis in the physical strip. However, 
at low temperature, this pattern will no longer provide the second largest 
eigenvalue. 

There is another definition of the transfer matrix for the i?-matrix (2.4), 
which is based on the graded formalism ]3|, ^, [TTJ, |2D[ of the quantum in- 
verse scattering method. In this case, the .R-matrix ( p.4|) is defined as 
R(v) = P g R(v), where P 9 is the graded permutation operator. As a re- 
sult, the ordinary row-to-row transfer matrix is defined as the super-trace of 
a monodromy matrix, and the phase factor (cf.[|16|, |T7|]) of the BAE disap- 
pears 0, p0| . In the QTM for the graded case, we are not sure how the phase 
factor corresponding to the right hand side in (|3.6| ) changes. 

One might suspect that our TBA equation (47?) is related to the one |63] 
for the Izergin-Korepin model since the representation space of osp( 1| 2) ^ 

(2) 

has close resemblance to that of the affine Lie algebra A 2 . However, to 
our knowledge, the TBA equation for the Izergin-Korepin model has been 
constructed |64| so that it reduces to the one for a sw(3)-invariant model in 
the rational limit. There will be another branch of the rational limit for the 
Izergin-Korepin model [ 17| . Together with the above-mentioned graded case, 



the comparison with our results will be interesting. 

It is an interesting problem to extend a similar analysis to more general 
osp{r\2s) model. Unfortunately, for r e Z> 2 case, we have Qonly the subset 
of the T-system [2B, While for r = 1 case, we have already a complete 
set of the T-system ||26|| , which will be relevant for the QTM method. So we 
hope to report, at the beginning, the TBA equation for the osp(l\2s) model 
in the near future. 

Finally we refer to another formulation of NLIE based on the finite num- 
ber of the auxiliary functions (see for example, 



53, 



5C| |7|). Al- 



though there seems to be no connection between the TBA equation and such 
NLIE, this formulation will be advantageous to the numerical calculation of 
the physical quantities at finite temperature. 



h osp(2\2) ~ sZ (1 1 2) case is an exception to this statement 



Thermodynamic Bethe ansatz equation 



16 



Acknowledgments 

ZT would like to thank Professor J. Suzuki for discussions at the earlier stage 



of the work pql , which gives a motivation of the present work. The authors 
also thank him for his useful comments. 

They would like to thank Professor A. Kuniba for his encouragement. 
This work is supported in part by a Grant-in- Aid for JSPS Research Fellows 
from the Ministry of Education, Science, Culture and Sports of Japan. 

References 

[1] V. Kac, Adv. Math. 26, 8 (1977). 

[2] P. P. Kulish and E. K. Sklyanin, J. Sov. Math 19, 1596 (1982). 

[3] P. P. Kulish, J. Sov. Math. 35, 2648 (1986). 

[4] V. V. Bazhanov and A. G. Shadrikov, Theor. Math. Phys. 73, 1302 
(1988). 

[5] P. P. Kulish and N. Yu. Reshetikhin, Lett. Math. Phys. 18, 143 (1989). 

[6] T. Deguchi, A. Fujii and K. Ito, Phys. Lett. B 238, 242 (1990). 

[7] H. Saleur, Nucl. Phys. B 336, 363 (1990). 

[8] R. B. Zhang, A. J. Bracken and M. D. Gould, Phys. Lett. B 257, 133 
(1991). 

[9] M. J. Martins and P. B. Ramos, J. Phys. A 27, L703 (1994). 

[10] P. Schlottmann, Phys. Rev. B 36, 5177 (1987). 

[11] F. H. L. Essler and V. E. Korepin, Phys. Rev. B 46, 9147 (1992). 

[12] F. H. L. Essler, V. E. Korepin and K. Schoutens, Phys. Rev. Lett. 68, 
2960 (1992). 

[13] A. Foerster and M. Karowski, Nucl. Phys. B 396, 611 (1993). 

[14] F. H. L. Essler and V. E. Korepin, Int. J. Mod. Phys. B 8, 3243 (1994). 

[15] Z. Maassarani, J. Phys. A 28, 1305 (1995). 

[16] M. J. Martins, Phys. Rev. Lett. 74, 3316 (1995). 



Thermodynamic Bethe ansatz equation 



17 



[17] M. J. Martins, Nucl. Phys. B 450, 768 (1995). 

[18] P. B. Ramos and M. J. Martins, Nucl. Phys. B 474, 678 (1996). 

[19] M. P. Pfannmuller and H. Frahm, Nucl. Phys. B 479, 575 (1996). 

[20] M. J. Martins and P. B. Ramos, Nucl. Phys. B 500, 579 (1997). 

[21] P. P. Kulish, N. Yu. Reshetikhin and E. K. Sklyanin Lett. Math. Phys. 
5, 393 (1981). 

[22] Z. Tsuboi, J. Phys. A 30, 7975 (1997). 

[23] Z. Tsuboi, Physica A 252, 565 (1998). 

[24] Z. Tsuboi, J. Phys. A 31, 5485 (1998). 

[25] Z. Tsuboi, Physica A 267, 173 (1999). 

[26] Z. Tsuboi, J. Phys. A 32, 7175 (1999). 

[27] N. Yu. Reshetikhin, Sov. Phys. JETP 57, 691 (1983). 

[28] N.Yu. Reshetikhin, Lett. Math. Phys. 14, 235 (1987). 

[29] V. V. Bazhanov and N.Yu. Reshetikhin, J. Phys. A 23, 1477 (1990). 

[30] J. Suzuki, Phys. Lett. A 195, 190 (1994). 

[31] A. Kuniba and J. Suzuki, Commun. Math. Phys. 173, 225 (1995). 

[32] A. Kuniba and J. Suzuki, J. Phys. A 28, 711 (1995). 

[33] A. Kuniba, Y. Ohta and J. Suzuki, J. Phys. A 28, 6211 (1995). 

[34] C.N. Yang and CP. Yang, J. Math. Phys. 10, 1115 (1969). 

[35] G. Jiittner, A. Kliimper and J. Suzuki, Nucl. Phys. B 512, 581 (1998). 

[36] H. Saleur, The continuum limit of sl(N/K) integrable super spin chains, 
Mv-int/9905007 . 

[37] M. J. Martins, B. Nienhuis and R. Rietman, Phys. Rev. Lett. 81, 504 
(1998). 

[38] K. Sakai and Z. Tsuboi, Thermodynamic Bethe ansatz equation for 
osp(l|2) integrable spin chain, [mat h-ph / 99 1 1 1 0| ; Mod. Phys. Lett. A, 
14, 2427 (1999). 



Thermodynamic Bethe ansatz equation 



18 



[39] M. Gaudin, Phys. Rev. Lett. 26, 1301 (1971). 

[40] M. Takahashi, Prog. Theor. Phys. 46, 401 (1971). 

[41] M. Takahashi and M. Suzuki, Prog. Theor. Phys. 48, 2187 (1972). 

[42] F.C. Alcaraz and M.J. Martins, J. Phys. A 21, L381, ibid. 4397 (1988). 

[43] H. L. Efiler, V. E. Korepin and K. Schoutens, J. Phys. A 25, 4115 (1992). 

[44] G. Jiittner and B.D. Dorfel, J. Phys. A 26, 3105 (1993). 

[45] M. Suzuki, Prog. Theor. Phys. 56, 1454 (1976). 

[46] M. Suzuki, Phys. Rev. B 31, 2957 (1985). 

[47] M. Suzuki and M. Inoue, Prog. Theor. Phys. 78, 787 (1987). 

[48] M. Inoue and M. Suzuki, Prog. Theor. Phys. 79, 645 (1988). 

[49] T. Koma, Prog. Theor. Phys. 78, 1213 (1987). 

[50] J. Suzuki, Y. Akutsu and M. Wadati, J. Phys. Soc. Jpn 59, 2667 (1990). 
[51] J. Suzuki, T. Nagao and M. Wadati, Int. J. Mod. Phys. B 6, 1119 (1992). 
[52] A. Kliimper, Ann. Physik 1, 540 (1992). 
[53] A. Kliimper, Z. Phys. B 91, 507 (1993). 

[54] G. Jiittner and A. Kliimper, Euro. Phys. Lett. 37, 335 (1997). 

[55] G. Jiittner, A. Kliimper and J. Suzuki, Nucl. Phys. B 487, 650 (1997). 

[56] A. Fujii and A. Kliimper, Nucl. Phys. B 546, 751 (1999). 

[57] K. Sakai, M. Shiroishi, J. Suzuki and Y. Umeno, Phys. Rev. B 60, 5186 
(1999). 

[58] A. Kuniba, K. Sakai and J. Suzuki, Nucl. Phys. B 525, 597 (1998). 
[59] J. Suzuki, Nucl. Phys. B 528, 683 (1998). 
[60] K. Sakai, J. Phys. Soc. Jpn. 68, 1789 (1999). 

[61] A. Kuniba, T. Nakanishi and J. Suzuki, Int. J. Mod. Phys. A 9, 5215 
(1994). 



Thermodynamic Bethe ansatz equation 



19 



[62] A. Kuniba, T. Nakanishi and J. Suzuki, Int. J. Mod. Phys. A 9, 5267 
(1994). 

[63] A. Kliimper and P. Pearce, Physica A 183, 304 (1992). 

[64] L. Mezincescu, R. I. Nepomechie, P. K. Townsend and A. M. Tsvelik, 
Nucl. Phys. B 406, 681 (1993). 



Thermodynamic Bethe ansatz equation 



List of Figures 





1 The distribution of the BAE roots which provides the largest 




eigenvalue for N = 12, u 


= 0.05. There are six two-string 


solutions which are symmetric with respect to the imaginary 


axis and the line = 3/4. 







2 Location of zeros for T^\(u,v) for m = 1,2,3, u = 0.05 and 


A = 12. Note that these zeros are distributed symmetrically 


with respect to the real and the imaginary axis 




3 Location of the roots of the BAE (3.6) for u = 0.05, N = 12 



and n — 11, which gives the second largest eigenvalue of the 
QTM. There are four two-strings and one three-string. They 
are symmetrically distributed with respect to the line = 
3/4 and the imaginary axis 





4 Location of zeros of t!£} 2 {u,v) for m — 1, 2,3, u — 0.05 and 


N = 12. Note that these zeros are distributed symmetrically 


with respect to the real and the imaginary axis. 






5 Location of the roots of the BAE (3.6) for u = 0.05, A = 14 



and n — 13, which gives the second largest eigenvalue of the 
QTM. There are six two-strings and one one-string. They are 
symmetrically distributed with respect to the line $sv = 3/4 
and the imaginary axis 





6 Location of zeros of T^ 9 (u,v) for m = 1,2, 3, u — 0.05 and 




A = 14. Note that these zeros are distributed symmetrically 






with respect to the real and the imaginary axis 



Thermodynamic Bethe ansatz equation 21 




Figure 1: The distribution of the BAE roots which provides the largest eigen- 
value for TV = 12, u — 0.05. There are six two-string solutions which are 
symmetric with respect to the imaginary axis and the line Ssv = 3/4. 
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Figure 2: Location of zeros for T^\(u,v) for m — 1,2,3, u — 0.05 and 
N = 12. Note that these zeros are distributed symmetrically with respect to 
the real and the imaginary axis. 
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Figure 3: Location of the roots of the BAE ( |3.6| ) for u = 0.05, N = 12 and 
n = 11, which gives the second largest eigenvalue of the QTM. There are 
four two-strings and one three-string. They are symmetrically distributed 
with respect to the line = 3/4 and the imaginary axis. 
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Figure 4: Location of zeros of T^\(u, v) for m — 1, 2, 3, u — 0.05 and iV = 12. 
Note that these zeros are distributed symmetrically with respect to the real 
and the imaginary axis. 
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Figure 5: Location of the roots of the BAE ( |3.6| ) for u = 0.05, N = 14 and 
n = 13, which gives the second largest eigenvalue of the QTM. There are 
six two-strings and one one-string. They are symmetrically distributed with 
respect to the line = 3/4 and the imaginary axis. 
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Figure 6: Location of zeros of T^\(u, v) for m — 1, 2, 3, u — 0.05 and iV = 14. 
Note that these zeros are distributed symmetrically with respect to the real 
and the imaginary axis. 
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